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ABSTRACT 

We have applied a new gauge-invariant, noncompact, Monte Carlo method 
to simulate U{1), SU (2), and SU{3) gauge theories on 8^ and 12^ lattices. The 
Creutz ratios of the Wilson loops agree with the exact results for U{1) for f3 > 1 
apart from a renormalization of the charge. The SU{2) and SU{3) Creutz 
ratios robustly display quark confinement at /3 = 0.5 and /3 = 1.5, respectively. 
At much weaker coupling, the SU{2) and SU{3) Creutz ratios agree with 
perturbation theory after a renormalization of the coupling constant. 

1. INTRODUCTION 

The first gauge-invariant noncompact simulations were carried out by Palumbo, 
Polikarpov, and Veselovlll and were based on earlier work by Palumbo et al.B They saw a 
confinement signal.^ Their action contains five terms, constructed from two invariants, 
and involves noncompact auxiliary fields and an adjustable parameter. 

The present paper implements and tests a new way of performing gauge-invariant 
noncompact simulations. This method is based upon a new noncompact action that 
is exactly invariant under lattice gauge transformations.^ The action is a natural dis- 
cretization of the classical Yang-Mills action with auxiliary fields that are compact 
group elements representing gauge transformations. 

We have used this method to simulate f/(l), SU{2), and SU{3) gauge theories 
on 8^ and 12^ lattices. The Creutz ratios of Wilson loops agree with the exact results 
for U{1) for (3 > 1 apart from a renormalization of the charge. The SU{2) and SU{3) 
Creutz ratios clearly show quark confinement at /3 = 0.5 and f3 = 1.5, respectively. 
At much weaker coupling, the SU (2) and SU (3) Creutz ratios agree with perturbation 
theory with a renormalized coupling constant. 

2. THE METHOD 

What constitutes a gauge transformation in this method? To find out, we look at 
the (massless) continuum fermionic action density iip'^ ^tf) . A suitable discretization 
of this quantity is iil){n)'^ ^[(l){n + e^) — il){n)\/a in which is a four-vector of integers 
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representing an arbitrary vertex of the lattice, is a unit vector in the /ith direction, 
and a is the lattice spacing. The product of Fermi fields at the same point is gauge 
invariant as it stands. The other product of Fermi fields becomes gauge invariant if we 
insert a matrix A^{n) of gauge fields 



-?/'(n)7^[(l + igaA^{n))i;{n + e^) - i){n)] 



(1) 



that transforms under a gauge transformation represented by the group elements U (n) 
and U{n + Cfj,) in such a way that 



1 + iagA'^{n) = U{n)[l + iagA^{n)]U'\n + e^). 
The required behavior is 



(2) 



A'^{n) = U{n)A^{n)U'\n + e^) + —U{n) \u-\n) - U-\n + e^) 



ag 

Let us define the lattice field strength F^y{n) as 

F^^{n) = ^[A^,{n + e^) - A^{n)] - ^[A^{n + e^) - A^{n)] 
+ ig[A^{n)A^{n + e^) - A^{n)A^{n + e^)] 

which reduces to the continuum Yang-Mills field strength in the limit a 
the gauge transformation (|^), this field strength transforms as 



(3) 



(4) 

0. Under 
(5) 



The field strength F^j_y{n) is antisymmetric in the indices /i and z/, but it is not hermitian. 
To make a positive plaquette action density, we use the Hilbert-Schmidt norm of Ffj,y{n) 



^ = ^Tr[Fj,(n)F,.(n)], 



(6) 



in which it is assumed that the generators of the gauge group have been orthonormalized 
as Tr(TaTfe) = k5ah- Because F^y{n) transforms covariantly (^, this action density is 
invariant under the noncompact gauge transformation (Q). 

The gauge transformation (|^) with group element U{n) = exp{—iaguj°'Ta) maps 
the usual matrix of gauge fields Af^{n) = Tay4^(n) into a matrix that generally lies 
outside the Lie algebra of the gauge group, although it does remain in the algebra to 
lowest (zeroth) order in the lattice spacing a. In our simulation we use for the gauge 
fields this more-general space of matrices. We use the action (^) in which the field 
strength (^) is defined in terms of gauge-field matrices that are of the form 



A^{n) = VAl{n)W-^ + —V (v-^ 

ag ^ 



w- 



(7) 



where A^^{n) is a matrix of gauge fields defined in the usual way, A°(n) = TaAl^^{n). 
Here the group elements V and W associated with the gauge field A^{n) are unrelated 
to those associated with the neighboring gauge fields A^{n + e^), A^{n), and A^{n + e^). 

We do not require the quantity l + igaA^{n) to be an element L^{n) of the gauge 
group. But if one did, then the matrix A^j_{n) of gauge fields would be related to the 
link L^{n) by A^{n) = {L^{n) — l)/{iga), and the action (|]) defined in terms of the 
field strength (HI) and this gauge-field matrix would be, mirabile dictu, Wilson's action: 



k - ^TiLf^{n)L„{n + e^)V^{n + eu)Ll[n 

RESULTS 



We have tested this method by applying it to the SU{2), and SU{?>) gauge 

theories on 8^ and 12^ lattices. In our initial configurations, the unitary matrices V 
and W were set equal to the identity matrix and the gauge fields A° were either zero 
or random. We allowed at least 10,000 sweeps for thermalization. 

For f/(l) and for /? > 1, our measured Creutz ratios agreed with the exact ones 
apart from finite-size effects and a renormalization of the charge. For instance at /3 = 1, 
we found on the 8^ lattice x(2, 2) = 0.142(1), x(2,3) = 0.098(1), x(3,3) = 0.047(1), 
X(2,4) = 0.085(1), x(3,4) = 0.030(1), and x(4,4) = 0.014(1). The first three of these 
x's are equal to the exact Creutz ratios for a renormalized value of I3r = 0.93; the last 
three are smaller than the exact ratios for (3^. = 0.93 due to finite-size effects by 6%, 
17%, and 42%, respectively. 

For SU{2) on the 8^ lattice at /? = 0.5, we found x(2,2) = 0.835(3), x(2, 3) = 
0.85(1), x(3, 3) = 0.9(2), and x(2,4) = 0.9(6) which robustly display confinement. At 
/? = 1, our six Creutz ratios track those of tree-level perturbation theory for a renor- 
malized value of = 1-75; the finite-size effects are hidden by incipient confinement. 

For SU{3) at (3 = 1.5, we found on the 8^ lattice x(2,2) = 1.175(3), x(2, 3) = 
1.16(2), x(3,3) = 1.7(16), and x(2,4) = 1.4(2), and on the 12^ lattice x(2, 2) = 
1.171(4), x(2,3) = 1.14(3), and x(2,4) = 1.9(7). At /3 = 2 we found on the 8^ lat- 
tice x(2,2) = 0.839(2), x(2,3) = 0.837(7), x(3,3) = 0.76(9), and x(2,4) = 0.87(3); 
and on the 12^ lattice x(2,2) = 0.832(2), x(2, 3) = 0.821(7), x(3, 3) = 0.71(7), and 
x(2,4) = 0.80(2). Within the limited statistics, these results exhibit confinement. At 
much weaker coupling, our ratios agree with perturbation theory apart from finite-size 
effects and after a renormalization of the coupling constant. 
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